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Abstract 

It is shown that the hadron prodnction in high energy pp, pp and e’''e“collisions, 
calcnlated by assnming that particles originate in hadron gas fireballs at thermal and 
partial chemical eqnilibrinm, agrees very well with the data. The temperature of 
the hadron gas hreballs, determined by htting hadron abundances, does not seem to 
depend on the centre of mass energy and kind of reaction, having a nearly constant 
value of about 170 MeV. This Ending supports a universal hadronization mechanism 
in all kinds of reactions consisting in a parton-hadron transition at critical values 
of temperature and pressure. 
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1 Introduction 


The thermodynamical approach to hadron production in hadronic collisions was originally 
introduced by Hagedorn |]^ almost thirty years ago. One of the most important phenomeno¬ 
logical indication of a thermal-like multihadron production in those reactions at high energy 
was found in the universal slope of transverse mass (i.e. rriT = -|- m?) spectra which 

is nowadays still being used extensively in heavy ions reactions, although it has been realized 
that tranverse collective motions of the hadron gas may signihcantly distort the primitive ther¬ 
mal m'T-spectrum [^. A much better probe of the existence of locally thermalized sources in 
hadronic collisions is the overall production rate of single hadron species which, being a Lorentz- 
invariant quantity, is not affected by local collective motions of the hadron gas. The analysis of 
hadron production rates with the thermodynamical ansatz implies that inter-species chemical 
equilibrium is attained, which is a tighter requirement with respect to a thermal-kinetical intra¬ 
species equilibrium assumed in the analysis of tut spectra. However, since hadron abundances 
can provide a cleaner indication of chemical and thermal equilibrium with respect to transverse 
mass spectra, we have focused our attention on the former issue. 

The smallness of the collision systems studied here requires appropriate theoretical tools: in 
order to properly compare theoretical predicted multiplicities to experimental ones, the use of 
statistical mechanics in its canonical form is mandatory, that means exact quantum numbers 
conservation is required, unlike in the grand-canonical formalism |Q. It will be shown indeed 
that particle average particle multiplicities in small systems are heavily affected by conservation 
laws well beyond what the use of chemical potentials predicts. 

The thermodynamical approach to hadron production in e’''e“collisions is here extended to 
pp and pp collisions; some assumptions which were made for e’''e“, in particular the restriction 
to two-jet events, are released; all calculations are performed with a bigger symmetry group 
(actually by also taking into account the conservation of electric charge). With this generalized 
model, we also update our results for e'''e“collisions. 

2 The model 

In ref. p, ^ a thermodynamical model of hadron production in e'''e“ collisions was developed 
on the basis of the following assumption: the hadronic jets observed in the hnal state of a 
e'''e“ —>■ qq event must be identihed with hadron gas phases having a collective motion. This 
identihcation is valid at the decoupling time, when hadrons stop interacting after their forma¬ 
tion and (possibly) a short expansion (freeze-out). Throughout this paper we will refer to such 
hadron gas phases with a collective motion as fireballs, following refs. |l|,P]. Since most events in 
a e’''e“ —>■ qq reaction are two-jet events, it was assumed that two hreballs are formed and that 
their internal properties, namely quantum numbers, are related to those of the corresponding 
primary quarks. In the so-called correlated jet scheme correlations between the quantum num¬ 
bers of the two hreballs were allowed beyond the simple correspondance between the hreball 
and the parent quark quantum numbers. This scheme turned out to be in better agreement 
with the data than a correlation-free scheme PJ. 

The more complicated structure of a hadronic collision does not allow a straightforward exten¬ 
sion of this model. If the assumption of hadron gas hreballs is maintained, the possibility of 
an arbitrary number of hreballs with an arbitrary conhguration of quantum numbers should 
be taken into account. To be specihc, let us dehne a vector Q = (Q, N, S, C, B) with integer 
components equal to the electric charge, baryon number, strangeness, charm and beauty respec- 
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tively. We assume that the hnal state of a pp or a pp interaction consists of a set of N hreballs, 
each with its own four-vector j3i = Ui/Ti, where Tj is the temperature and u, = ( 7 *, Ay*) is the 
four-velocity 0], quantum numbers Q° and volume in the rest frame V^. The quantum vectors 
Q° must fulhll the overall conservation constraint = Q° where Q° is the vector of the 

initial quantum numbers, that is Q° = ( 2 , 2 , 0 , 0 , 0 ) in a pp collision and Q° = ( 0 , 0 , 0 , 0 , 0 ) in 
a pp collision. 

The invariant partition function of a single hreball is, by dehnition: 


z.(Ql) = E 

states 


3 


Qi,Q° 


( 1 ) 


where Pi is its total four-momentum. The factor ^q.,qo is the usual Kronecker tensor, which 
forces the sum to be performed only over the hreball states whose quantum numbers Q* are 
equal to the particular set Q°. It is worth emphasizing that this partition function corresponds 
to the canonical ensemble of statistical mechanics since only the states fulhlling a hxed chemical 
requirement, as expressed by the factor ^q.^qo, are involved in the sum ( 1 ). 

By using the integral representation of Sq. qo: 


27r27r27r27r27r 




(27r) 


0 Jo Jo Jo Jo 




Eq. (1) becomes: 


( 2 ) 


= E 

states 


1 






( 3 ) 


This equation could also have been derived from the general expression of partition function 
of systems with internal symmetry P, ^ by requiring a U(l)® symmetry group, each U(l) 
corresponding to a conserved quantum number; that was the procedure taken in ref. |p. 

The sum over states in Eq. 3 can be worked out quite straightforwardly for a hadron gas of 
Nb boson species and Np fermion species. A state is specihed by a set of occupation numbers 
{uj^k} for each phase space cell k and for each particle species j. Since Pi = J2j,kPknj,k and 
Qi = J2j,k^j'^j,k, where = {Qj, Nj, Sj,Cj, Bj) is the quantum number vector associated to 
the particle species, the partition function (3) reads, after summing over states: 


1 J'B 1\p 

Z^iQ'^) = 77 ^ [ dVe‘^^^exp ^ log (1 - ^ log (1 + . 

j=i k j=i k 

( 4 ) 

The last expression of the partition function is manifestly Lorentz-invariant because the sum 
over phase space is a Lorentz-invariant operation which can be performed in any frame. The 
most suitable one is the hreball rest frame, where the four-vector jSi reduces to: 

A = (;^, 0,0,0) (5) 

Ti being the temperature of the hreball. Moreover, the sum over phase space cells in Eq. (4) 
can be turned into an integration over momentum space going to the continuum limit: 
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where V is the hreball volume and Jj the spin of the hadron. As in previous studies on 
e'''e“collisions and heavy ions collisions [|T^], we supplement the ordinary statistical me¬ 
chanics formalism with a strangeness suppression factor 7 ^ accounting for a partial strangeness 
phase space saturation; actually the Boltzmann factor of any hadron species containing 

s strange valence quarks or antiquarks is multiplied by yf. With the transformation ( 6 ) and 
choosing the hreball rest frame to perform the integration, the sum over phase space in Eq. (4) 
becomes: 


log (1 ± 


(2 J, + 1) E, 


(27r) 


- J d^p log (1 ± ^ V, Fj{T„ 7 ,) 


(7) 


where the upper sign is for fermions, the lower for bosons and Vi is the hreball volume in its 
rest frame; the function Fj{Ti, 7 ^) is a shorthand notation of the momentum integral in Eq. (7). 
Hence, the partition function (4) can be written: 


Z.(Q?)= d',#.e>«?-<'expmx:nn..7.)l. (8) 

The mean number {nj)i of the particle species in the hreball can be derived from 2’(Q?) 
by multiplying the Boltzmann factor exp {—^Jp^ -|- mj/T), in the function Fj in Eq. ( 8 ) by a 
hctitious fugacity \j and taking the derivative of logZj(Q°, Xj) with respect to Xj at Xj = 1 : 


{nj)i = ^logZi(Q°,Aj) 


A,=l 


(9) 


The partition function Zj(Qj,Aj) supplemented with the Xj factor is still a Lorentz-invariant 
quantity and so is the mean number (n^)*. From a more physical point of view, this means that 
the average multiplicity of any hadron does not depend on hreball collective motion, unlike its 
mean number in a particular momentum state. 

The overall average multiplicity of the hadron, for a set of N hreballs in a certain quantum 
conhguration {Q?,..., 0,%} is the sum of all mean numbers of that hadron in each hreball: 


N 


d 


K) = 


i=l 


d 


N 


iogn^.(Q°,A,) 


2=1 


A,=l 


( 10 ) 


In general, as the quantum number conhgurations may huctuate, hadron production should be 
further averaged over all possible hreballs conhgurations Q ^,fulhlling the constraint 
Eh Q° = Q°. To this end, suitable weights w(Q?,..., Qn), representing the probability of 
conhguration {Q^,..., Q^} to occur for a set of N hreballs, must be introduced. Basic features 
of those weights are: 


N 

w(Q°,...,Q%)=0 it 

2=1 

«>(q;,...,q 7 = i. ( 11 ) 

q!.qS, 

For the overall average multiplicity of hadron j we get: 
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( 12 ) 


((%)) = 


E 

q?,...,qS 




d 

d\j 


N 


iogn^.(Q°,A,; 


2=1 


There are infinitely many possible choices of the weights w{Q^,, Q%), all of them equally 
legitimate. However, one of them is the most pertinent from the statistical mechanics point of 
view, namely: 


w(q;.---.Qk) = 


Se.qiq<. n.., zm 


(13) 


I]Q0,...,Q0^SiQ9,Q0 Hill ) 

It can be shown indeed that this choice corresponds to the minimal deviation from statistical 
equilibrium of the system as a whole. In fact, putting weights (13) in the Eq. (12), one obtains: 


d 


N 


%)) = E '^S,Q0,Q0 Aj- 


2=1 


A,=l 


(14) 


This means that the average multiplicity of any hadron can be derived from the following 
function of Q°: 


^(Q°) 


= E^ 


SiQ°,QO 


Q?,.,Q 


0 

N 


n z.(Q?). 


2=1 


(15) 


with the same recipe given for a single fireball in Eq. (9). By using expression (1) for the 
partition functions Zj(Q°), Eq. (15) becomes: 


Since 


N 

^(Q°) = E n E 

Q5,...,Q0, i=l states; 



(16) 


the function (16) can be written as 

Z(Q“) =E ■ ■ ■ E . (18) 

statesi statesN 

This expression demonstrates that Z{Q^) may be properly called the global partition function 
of a system split into N subsystems which are in mutual chemical equilibrium but not in mutual 
thermal and mechanical equilibrium. Indeed it is a Lorentz-invariant quantity and, in case of 
complete equilibrium, i.e. Pi = P2 = ■ ■ ■ = Pn = P, it would reduce to: 

Z(Q“) = E ... E -'’»>5aQ.,Q. = E e-« '’5q,Q. . (19) 

states 1 states^ states 

which is the basic definition of the partition function. 

To summarize, the choice of weights (13) allows the construction of a system which is out of 
equilibrium only by virtue of its subdivision into several parts having different temperatures 
and velocities. Another very important consequence of that choice is the following: if we assume 
that the freeze-out temperature of the various fireballs is constant, that is Ti = ... = T^r = T, 
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and that the strangeness snppression factor % is constant too, then the global partition fnnction 
(18) has the following expression: 


Z(Q°) = / d".#. e ■ exp |(E.l/.) F, (T. ^,)] . (20) 

Here the V)’s are the hreball volnmes in their own rest frames; a proof of (20) is given in 
ref. 1^. Eq. (20) demonstrates that the global partition fnnction has the same fnnctional form 
(3), (4), ( 8 ) as the partition fnnction of a single hreball, once the volnme Vi is replaced by the 
global volume V = Vi. Note that the global volnme absorbs any dependence of the global 
partition fnnction (20) on the nnmber of hreballs N. Thns, possible variations of the nnmber 
N and the size Vi of hreballs on an event by event basis can be tnrned into hnctuations of 
the global volnme. In the remainder of this section these hnctnations will be ignored; in Sect. 
3 it will be shown that they do not ahect any of the following resnlts on the average hadron 
mnltiplicities. 

The average mnltiplicity of the hadron can be determined with the formulae (14)-(15), by 
using expression (20) for the function Z(Q°): 

y d^ e ‘ exp [H ^ Fj (T, 7 ,)] 

X ^ 

(27r)3 J gxp (y/p^ + rrij/T + icp- • 0 ) ± 1 

where the upper sign is for fermions and the lower for bosons. This formula can be written in 
a more compact form as a series: 


W /'-.-I ^^l+l ... ■2^(Q° ^Hj) 

m)) = Z^^(Tl) Is - , (22) 

where the functions Zj{ri) are dehned as: 

Zj^n) = {2Jj + 1 ) J d^p exp + m]/T) = {2Jj + 1 ) K 2 (^) . (23) 

K 2 is the McDonald function of order 2. Eq. (22) is the hnal expression for the average 
multiplicity of hadrons at freeze-out. Accordingly, the production rate of a hadron species 
depends only on its spin, mass, quantum numbers and strange quark content. The chemical 
factors Z(Q°—nqj )/Z(Q°) in Eq. (22) are a typical feature of the canonical approach due to the 
requirement of exact conservation of the initial set of quantum numbers. These factors suppress 
or enhance production of particles according to the vicinity of their quantum numbers to the 
initial Q° vector. The behaviour of Z{Cl) as a function of electric charge, baryon number and 
strangeness for suitable T, V and 7 ^ values is shown in Fig. 1; for instance, it is evident that 
the baryon chemical factors Z( 0 , A^, 0 , 0 , 0 )/Z( 0 , 0 , 0 , 0 , 0 ) connected with an initially neutral 
system play a major role in determining the baryon multiplicities. The ultimate physical 
reason of “charged” particle (q^ 7 ^ 0 ) suppression with respect to “neutral” ones (q^ = 0 ), in 
a completely neutral system (Q° = 0), is the necessity, once a “charged” particle is created, 
of a simultaneous creation of an anti-charged particle in order to fulhll the conservation laws. 
In a finite system this pair creation mechanism is the more unlikely the more massive is the 


6 












lightest particle needed to compensate the hrst particle’s quantum numbers. For instance, once 
a baryon is created, at least one anti-nucleon must be generated, which is rather unlikely since 
its mass is much greater than the temperature and the total energy is hnite. On the other hand, 
if a non-strange charged meson is generated, just a pion is needed to balance the total electric 
charge; its creation is clearly a less unlikely event with respect to the creation of a baryon as 
the energy to be spent is lower. This argument illustrates why the dependence of 2’(Q) on 
the electric charge is much milder that on baryon number and strangeness (see Fig. 1). In 
view of that, the dependence of Z{Q) on electric charge was neglected in the previous study on 
hadron production in e’''e“collisions [^. These chemical suppression effects are not accountable 
in a grand-canonical framework; in fact, in a completely neutral system, all chemical potentials 
should be set to zero and consequently “charged” particles do not undergo any suppression 
with respect to “neutral” ones. 


3 Other constraints and approximations 


Actually, the global partition function (18) has to be modihed in pp collisions owing to a major 
effect in such reactions, the leading baryon effect |TT[. Indeed, the sum (18) includes states with 
vanishing net absolute value of baryon number, whereas in pp collisions at least one baryon- 
antibaryon pair is always observed. Hence, the simplest way to account for the leading baryon 
effect is to exclude those states from the sum. Thus, if |iV| = 1]^ |iVj| denotes the absolute value 
of the baryon number of the system, the global partition function (18) should be turned into: 


^=E 


E e 

states 1 statesfvf 


-dl-Pl Q~dN-PN^. 


SiQi.QO 


■ E--- E 

statesi statesfvf 




(24) 


The hrst term, that we dehne as ^i(Q°), is equal to the function Z{Qf) in Eqs. (18), (20), 
while the second term is the sum over all states having vanishing net absolute value of baryon 
number. The absolute value of baryon number can be treated as a new independent quantum 
number so that the processing of the partition function described in Eqs. (l)-(3) can be repeated 
for the second term in Eq. (24) with a U(l)® symmetry group. Accordingly, this term can be 
naturally denoted by Z2(Q°,0), so that Eq. (24) reads: 


Z = Zi(Q°)-Z2(Q°,0) . (25) 

By using the integral representation of (5|Ar|,o 

Vl,o = 7^/ dV^e’l^l-^ (26) 

2n Jo 

in the second term of Eq. (24), one gets: 


Z2(Q",0) = ^y'dye‘«“-*explV'i:F,(r,7.)] 

X J di- expE 1 dy log (1 + , (27) 

where the hrst sum over j runs over all mesons and the second over all baryons. The average 
multiplicity of any hadron species can be derived from the global partition function (25) with 
the usual prescription: 
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(28) 




The calculation of the average multiplicity of primaries according to Eq. (22) (and Eq. (28)) 
involves several rather complicated Eve-dimensional integrals which have been calculated nu¬ 
merically after some useful approximations, described in the following. Since the temperature 
is expected to be below 200 MeV, the primary production rate of all hadrons, except pions, is 
very well approximated by the first term of the series ( 22 ): 


({%» - 7 ? Zj 


Z(ce - q,) 

Z(Q») 


(29) 


where we have put Zj = Zj{i). This approximation corresponds to the Boltzmann limit of Fermi 
and Bose statistics. Actually, for a temperature of 170 MeV, the primary production rate of 
K"*", the lightest hadron after pions, differs at most (i.e. without the strangeness suppression 
parameter and the chemical factors which further reduce the contribution of neglected terms) 
by 1.5% from that calculated with Eq. (29), well within usual experimental uncertainties. 
Corresponding Boltzmannian approximations can be made in the function Z{Q), namely 


log (1 ± ~ ^ 

which turns Eq. (20) (for a generic Q) into: 


z(q)~ 


(2i) 


d^0 


PQ- 


exp 


E 






+ E 


V 


(27r)= 


d^p log (1 - ^ 

(31) 


where the first sum runs over all hadrons except pions and the second over pions. 

As a further consequence of the expected temperature value, the functions of all charmed 
and bottomed hadrons are very small: with T = 170 MeV and a primary production rate of K 
mesons of the order of one, as the data states, the function of the lightest charmed hadron, 
D°, turns out to be ~ 10“^; chemical factors produce a further suppression of a factor 10“^. 
Therefore, thermal production of heavy flavoured hadrons can be neglected, as well as their 
functions in the exponentiated sum in Eq. (31), so that the integration over the variables 04 
and 05 can be performed: 


Z{Q,C,B) ~ ^|dye'«-«expE2,77e-‘'l’'« 

- E 7%3 /d’P log (1 - ,32) 

j = l 

Q and are now three-dimensional vectors consisting of electric charge, baryon number, and 
strangeness; the five-dimensional integrals have been reduced to three-dimensional ones. 

Apart from the hadronization contribution, which is expected to be negligible in this model, 
production of heavy flavoured hadrons in high energy collisions mainly proceeds from hard 
perturbative QCD processes of cc and bb pairs creation. The fact that promptly generated 
heavy quarks do not reannihilate into light quarks indicates a strong deviation from statistical 
equilibrium of charm and beauty, much stronger than the strangeness suppression linked with 



Nevertheless, it has been found in e+e“collisions that the relative abundances of charmed 
and bottomed hadrons are in agreement with those predicted by the statistical equilibrium 
assumption, conhrming its full validity for light quarks and quantum numbers associated to 
them. The additional source of heavy flavoured hadrons arising from perturbative processes 
can be accounted for by modifying the partition function (31). In particular, the presence of 
one heavy flavoured hadron and one anti-flavoured hadron should be demanded in a fraction 
of events = cr(/ ^ cc)/(j(I) (or i?;, = a{I —> bb)/cr(I)) where / is the initial colliding 
system, namely e’''e“, pp or pp , and cr(J) is meant to be the hadronic, inelastic and non-single- 
diffractive cross section respectively. Accordingly, the partition function to be used in events 
with a perturbative cc pair, is, by analogy with Eq. (24)-(25) and the leading baryon effect: 

^ = E---E 

statesi states^ 

- ■■■e”'’"‘'’'''5s.q.,Q««|C|.o = Zi(Q”)-Z2(Q",0), (33) 

statesi statesN 

where ICI is the absolute value of charm. The partition function to be used in events with a 
perturbative bb is, mutatis mutandis, analogous. Therefore, in the following analysis, we limit 
ourselfs to the case of a cc pair, being understood that all derived results hold for the bb case 
as well. 

In order to derive the average multiplicity of charmed hadrons in events in which one cc pair 
is created owing to a hard QCD process, the usual Eq. (28) should be used with the partition 
function Z in Eq. (33). The Z can be worked out the same way as for the leading baryon effect, 
namely: 


Z = Zi(QO)-Z2(Q°,0). 


where the function Zi can be written as in Eq. (31): 


(34) 


^i(Q) ^ 


1 

(27r)® 


j 



V 

(27r)3 


and the function Z 2 : 


(35) 


Z 2 (Q“.A') = 

+ E 7X^5/dy log (36) 

j=i (27r)3 J 

where the second sum in the exponentials in both Eqs. (35) and (36) runs over the three pion 
states and \Cj\ in Eq. (36) is the absolute value of hadron’s charm. 

Henceforth Q° and q*, q^, q^ are to be understood as three-dimensional vectors having as 
components electric charge, baryon number and strangeness, while charm and beauty will be 
explicitely written down. By using this notation, the average multiplicity of a charmed hadron 
with Cj = 1 turns out to be (see Eq. (29)): 


9 








(37) 


// w _ ^i(Q° - qj',-1,0) - ^2(Q°0,-1) 

Zi(QO,0,0)-Z2(QO, 0,0,0) 

Since the 2 : functions of heavy flavoured hadrons are 1 a power expansion in the Zj'’s of all 
charmed and anticharmed hadrons can be performed from Zj = 0 in the integrands of Eqs. (35) 
and (36), that is: 


expE 


Is ^3^' 


-ICiyi 


~ 1 + 


E 

3 




-iqj-0 I ± 


+ 2!^^ s'Ys^^iZje 

i,3 




for Eq. (35) and 


(38) 


exp 


E 

3 




Ih] ^ 


1 + E ^3^' 


^ i,3 


s 'Is^ZiZje 




(39) 


for Eq. (36). Furthermore, the functions of the bottomed hadrons can be negelected as they 
are 1 as well and beauty in Eq. (37) is always set to zero. 

Those expansions permit carrying out integrations in the variables ip, 04 and 05 in Eqs. (35) 
and (36). Thus: 


^i(Q° - qj, -1,0) ^ Y.^s'ZiC{Cf - qj- - qi) , (40) 

i 

where the sum runs over the anticharmed hadrons as the integration in 04 of terms associated to 
charmed hadrons yields zero. The ( function on the right-hand side is the same as in Eq. (32). 
Moreover: 


z,(Q», 0,0) ~ C(Q") + E 7?7'*2.2»C(Q” - q* - qi), (4i) 

i,k 

where the index i runs over all charmed hadrons and index k over all anti-charmed hadrons. 
Owing to the presence of the absolute value of charm in the exponential exp[ 110^10] in its 
integrand function, the function Z 2 {QP, C, B, K) vanishes if it' < 0 and yields K^^-oidei terms 
of the power expansion in Zj if iT > 0 (see Eq. (36)). Therefore: 


Z2(q°,o,o,o) = C(Q°) 

(42) 

and 


Z2(Q°,-1,0,-1) = 0 . 

(43) 

Finally, inserting Eqs. (40), (41), (42) and (43) in Eq. (37) one gets: 


^ ^ Ei,fc7078"2:i^fcC(Q°-qi-qfc) ’ 

(44) 


where the indices j, k label charmed hadrons and i labels anticharmed hadrons. From previous 
equation it results that the overall number of charmed hadrons is 1 , as it must be if c quark 
production from fragmentation is negligible. The average multiplicity of anticharmed hadrons 
is of course equal to charmed hadrons one. The same formula (44) holds for the average 
multiplicity of bottomed hadrons in events with a perturbatively generated bb pair. If leading 
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baryon effect is taken into account for pp collisions, the formula (44) gets more complicated, 
but the procedure is essentially the same. 

It is clear that a possible charm or beauty suppression parameter 7 c or introduced by 
analogy with strangeness suppression parameter 7 ^, would not be revealed from the study of 
heavy flavoured hadron production because a single factor multiplying all Zj functions would 
cancel from the ratio in the right-hand side of Eq. (44). 

So far, we tacitly assumed that the parameters T, V and 7 ^ do not fluctuate on an event by 
event basis. If freeze-out occurs at a hxed hadronic density in all events, then it is a reasonable 
ansatz that T and 7 ^ do not undergo any fluctuation since the density mainly depends on those 
two variables. However, there could still be signihcant volume fluctuations due to event by 
event variations of the number and size of the hreballs from which the primary hadrons emerge 
(see end of Sect. 2). However, it can be shown that, as far as the average hadron multiplicity 
ratios are concerned, possible fluctuations of V are negligible. To this end, let us dehne piV) 
as the probability density of picking a volume between V and V + dV in a single event. The 
primary average multiplicity of the hadron is then; 

{(n,» = j AVp{V) |;(Tl)”+‘ ■ (45) 

Provided that the volume V does not fluctuate over a too large range, the chemical factors can 
be taken out of the integral in Eq. (45) and evaluated at the mean volume V, because their 
dependence on it is mild and can be neglected (see Fig. 2). Therefore, the integrand depends 
on the volume only through the functions Zj(^n) whose dependence on V is linear (see Eq. (23)) 
and which can be reexpressed as 


^j{n){V, T, 7 ,) = V ^j(n)(T, %) . (46) 

Then, from Eq. (45), 

(WA) ^ E(=fi)"+‘ 7r'?j(«)(r. 7 .)/dvxr)r. (47) 

The integral on the right-hand side is the mean volume V. Thus: 

{{«,)) f;(Tl)’‘+‘7r'.-yw(V.r,7,) . (48) 

It turns out that the relative hadron abundances do not depend on the volume fluctuations 
since the mean volume V appearing in the above equation (replacing the volume V in Eq. (23)) 
is the same for all species: all results previously obtained are unaffected. 

4 Fit procedure and results 

The model described so far has three free parameters: the temperature T, the global volume 
V and the strangeness suppression parameter 7 ^. They have to be determined by a £t to the 
available data on hadron inclusive production at each centre of mass energy. Eq. (28) yields 
the mean number of hadrons emerging directly from the thermal source at freeze-out, the so- 
called primary hadrons, iHl as a function of the three free parameters. After freeze-out, 
primary hadrons trigger a decay chain process which must be properly taken into account 
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in a comparison between model predictions and experimental data, as the latter generally 
embodies both primary hadrons and hadrons generated by heavier particles decays. Therefore, 
in order to calculate overall average multiplicities to be compared with experimental data, the 
primary yield of each hadron species is added to the contribution stemming from the decay of 
heavier hadrons, which is calculated by using experimentally known decay modes and branching 
ratios 
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The primary average multiplicity is calculated by using Eq. (22) for pp collisions and Eq. (28) 
for pp collisions. For e^e“collisions the primary production rate of the hadron is given by: 


d 


{(nj» = i:B,g-logZ. 


i=l 


(49) 


where: 


Ri = 


cr e^e" 




(Ti^e'^e —hadrons) ^ ^ 

are the branching fractions into specihc flavours and Zj are the corresponding partition func¬ 
tions, namely Eq. (32) for u, d quarks, Eq. (33) for c quarks and also for s and b quarks where 
\C\ is replaced by |S| and [S'! respectively. The fractions Ri depend on the centre of mass en¬ 
ergy of the colliding system and determine the primary production of heavy flavoured hadrons 
to be taken into account as explained in the previous section. 

As far as hadronic collisions are concerned, the fractions Rc and Rb are worst known with 


respect to e“''e collisions. Available data on charm cross-sections indicate a fraction 


Rc ~ 10“^ 10“^ at centre of mass energies < 30 GeV and, consequently, much lower val¬ 

ues for bottom quark production. Therefore, the perturbative production of heavy quarks in 
hadronic collisions can be neglected as long as one deals with light flavoured hadron production 
at a/s < 30 GeV. We assume that it may be neglected at any centre of mass energy. An upper 
estimate of charm production in pp collisions at 900 GeV {Rc = 0.3) based on theoretical calcu¬ 


lations [0, does not affect signihcantly the htted parameters. All light flavoured hadrons and 
resonances with a mass < 1.7 GeV have been included among the primary generated hadron 
species. All heavy flavoured hadrons included in the JETSET 7.4 [|^ program tables have 
been used. 

The mass of resonances with T > 1 MeV has been distributed according to a relativistic Breit- 
Wigner function within ±2r from the central value. The y* strangeness suppression factor has 
also been applied to neutral mesons such as /, a;, etc. according to the their strange valence 
quark content; mixing angles quoted in [13| have been used. Once the average multiplicities 


of the primary hadrons have been calculated as a function of the three parameters T, V and 
js, the decay chain is performed until vr, /i, K^, K°, A, S, 0“ or stable particles are 
reached, in order to match the average multiplicity dehnition in pp and pp collisions experi¬ 
ments. The decay chain is further extended in e’''e“collisions since experiments also include 
the decay products of K^, A, S, and 0“ in their multiplicity dehnition. Finally, the overall 
yield is compared with experimental measurements, and the 




= ^(theoj — expe/^/error^ 


(51) 


is minimized. 

As far as the data set is concerned, we used all available measurements of hadron multiplicities 
in non-single-diffractive pp [|^ and inelastic pp collisions down to a centre of mass energy 
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of about 19 GeV and e’''e collisions between 29 and 91.2 GeV 0. 


Whenever several measurements at the same centre of mass energy have been available, aver¬ 
ages have been calculated according to a weighting procedure described in ref. prescribing 
rescaling of errors to take into account a posteriori correlations and disagreements of experi¬ 
mental results. If an experiment repeated the same measurement at the same centre of mass 
energy, we used only the most recent value. 

Since the decay chain is an essential step of the htting procedure, calculated theoretical multi¬ 
plicities are affected by experimental uncertainties on masses, widths and branching ratios of 
all involved hadron species. In order to estimate the effect of these uncertainties on the results 
of the fit, a two-step procedure for the fit itself has been adopted: firstly, the fit has been 
performed with a including only experimental errors and a set of parameters Tq, Vq, has 
been obtained. Then, the various masses, widths and branching ratios have been varied in turn 
by their errors, as quoted in ref. |T^, and new theoretical multiplicities calculated, keeping the 
parameters Tq, Vq, 7so fixed. The differences between old and new theoretical multiplicity values 
have been considered as additional systematic errors to be added in quadrature to experimental 
errors. Finally, the fit has been repeated with a including overall errors so as to obtain final 
values for model parameters and for theoretical multiplicities. Among the mass, width and 
branching ratio uncertainties, only those producing significant variations of final hadron yields 
(actually more than 130) have been considered. 

The fitted values of the parameters T, V, % at various centre of mass energy points are quoted 
in Table 1. The fit quality is very good at almost all centre of mass energies as demonstrated 
by the low values of y^’s and by the Figs. 3, 4, 5, 6, 7. Owing to the relatively large value of 
at ^/s = 27.4 GeV in pp collisions and at ^/s = 91.2 GeV in e’''e“collisions, variations of 
fitted parameters larger than fit errors must be expected when repeating the fit excluding data 
points with the largest deviations from the theoretical values. Therefore, the fit at ^/s = 27.4 
GeV pp collisions has been repeated excluding in turn (A°, 0) and (K“, pions), respectively, 

from the data set; the fit at = 91.2 GeV e’''e“collisions has been repeated excluding E*^. 
The maximum difference between the new and old fit parameters has been considered as an 
additional systematic error and is quoted in Table 1 within brackets. The fitted temperatures 
are compatible with a constant value at freeze-out independently of collision energy and kind 
of reaction (see Fig. 8) and it is in good agreement with that found in heavy ions collisions 
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The 7 s parameter exhibits a very slow rise from 20 to 900 GeV (see Fig. 9); its value of ~ 0.5 
over the whole explored centre of mass energy range in hadronic collisions proves that complete 
strangeness equilibrium is not attained. Furthermore, 7 s appears to be definitely lower in pp 
and pp collisions than in e''“e“collisions at equal centre of mass energy. 

On the other hand, the global volume does increase as a function of centre of mass energy as 
it is proportional, for nearly constant T and 7 s, to overall multiplicity which indeed increases 
with energy. Its values range from 6.4 fm^ at a/s = 19.4 GeV pp collisions, at a temperature 
of 191 MeV, up to 67 fm^ at y/s = 900 GeV pp collisions at a temperature of 170 MeV. 

Once T, V and jg are determined by fitting average multiplicities of some hadron species, their 
values can be used to predict average multiplicities of any other species, at a given centre of 
mass energy. 

The physical significance of the results found so far depends on their stability as a function 
of the various approximations and assumptions which have been introduced. First, the tem¬ 
perature and 7 s values are low enough to justify the use of the Boltzmann limits (29), (30) 
for all hadrons except pions, as explained in Sect. 3. As far as the effect of a cut-off in the 
hadronic mass spectrum goes, the most relevant test proving that our results do not depend 
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Table 1: Values of fitted parameters in pp, pp and e’''e collisions. 


(GeV) 

T(MeV) 

VT3 

7. 

xVdof 

pp collisions 

19.5 

190.8 ±27.4 

5.8 ±3.1 

0.463 ±0.037 

6.38/4 

23.8 

194.4 ± 17.3 

6.3 ±2.5 

0.460 ±0.067 

2.43/2 

26.0 

159.0 ±9.5 

13.4 ±2.7 

0.570 ±0.030 

1 .86/2 

27.5 

169.0 ±2.1 

11.04 ±0.69 

0.510 ±0.011 

136.4/27 


(±3.4) 

(±1.4) 

(±0.025) 



pp collisions 


200 

175.4 ± 14.8 

24.3 ± 7.9 

0.537 ±0.066 

0.698/2 

546 

181.7 ± 17.7 

28.5 ± 10.4 

0.557 ±0.052 

3.80/1 

900 

170.2 ± 11.8 

43.2 ± 11.8 

0.578 ±0.063 

1.79/2 


e’''e collisions 


29^-30 

163.6 ±3.6 

15.2 ± 1.4 

0.724 ±0.045 

24.7/13 

34^-35 

165.2 ±4.4 

14.6 ± 1.6 

0.788 ± 0.045 

10.5/8 

42.6 H- 44 

169.6 ±9.5 

14.7 ±3.1 

0.730 ±0.060 

4.94/4 

91.2 

160.6 ± 1.7 

26.4 ± 1.3 

0.675 ±0.020 

60.8/21 


(±3.1) 

(±1.1) 

(±0.027) 



on it is the stability of the number of primary hadrons against changes of the cut-off mass. 
The £t procedure intrinsically attempts to reproduce fixed experimental multiplicities; if the 
number of primary hadrons does not change significantly by repeating the fit with a slightly 
lower cut-off, the production of heavier hadrons excluded by the cut-off must be negligible, in 
particular with regard to its decay contributions to light hadron yields. In this spirit all fits 
have been repeated moving the mass cut-off value from 1.7 down to 1.3 GeV in steps of 0.1 
GeV, checking the stability of the amount of primary hadrons as well as of the fit parameters. 
It is worth remarking that the number of hadronic states with a mass between 1.7 and 1.6 GeV 
is 238 out of 535 overall, so that their exclusion is really a severe test for the reliability of the 
final results. Figure 10 shows the model parameters and the primary hadrons in pp collisions 
at y/s = 900 GeV; above a cut-off of 1.5 GeV the number of primary hadrons settles at an 
asymptotically stable value, whilst the fitted values for T, V, % do not show any particular 
dependence on the cut-off. Therefore, we conclude that the chosen value of 1.7 GeV ensures 
that the obtained results are meaningful. 

5 Conclusions 

A thermodynamical approach to hadronization, based on the use of the canonical formalism of 
statistical mechanics ensuring the exact quantum numbers conservation, succeeds in reproduc¬ 
ing the average multiplicities of hadron species in e’''e“, pp and pp collisions. These quantities 
are best suited in investigating the main characteristics of hadronization because of their in¬ 
dependence from local collective flows. The success of this statistical approach indicates that. 
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apart from the violation of strangeness phase space saturation as implemented in the 7 ^ pa¬ 
rameter, hadron production in elementary high energy collisions is dominated by phase space 
rather than by dynamical effects. The dynamics is essentially contained in hadron gas hreballs 
collective motions reflecting local hard partons kinematics. 

The temperature at freeze-out seems to be constant independently of kind of reaction and cen¬ 
tre of mass energy. This uniformity suggests that parton-hadron transition cannot occur before 
the parameters of prehadronic matter such as density or pressure have dropped below critical 
values corresponding to a temperature of about 170 MeV in a (partially) equilibrated hadron 
gas. This value is in agreement with lattice QCD results . 


From previous results it turns out that strangeness suppression, as well as the survival of per- 
turbatively created cc and bb pairs are the only trace to strong interaction dynamics before 
hadronization. This hnding suggests that such non-equilibrium effects are mainly related to 
quark mass thresholds. 
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Figure 5 
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Figure 7 


Figure 8 


Figure 9 


Figure 10 


Behaviour of the global partition function Z as a function of electric charge, baryon 
number and strangeness, keeping all remaining quantum numbers set to zero, for T = 170 
MeV, V = 20 fm^ and 7 ^ = 0.5 

Behaviour of the non-strange baryon chemical factor Z(0,1, 0, 0, 0)/Z(0, 0, 0, 0, 0) as a 
function of volume for different values of the temperature and a hxed value 7 ^ = 0.5 for 
the strangeness suppression parameter 

Results of hadron multiplicity hts for pp collisions at centre of mass energies 19.4, 23.8 
and 26 GeV. Experimental average multiplicities are plotted versus calculated ones. The 
dashed lines are the quadrant bisectors: well htted points tend to lie on these lines 

Results of hadron multiplicity £t for pp collisions at centre of mass energy of 27.4 GeV. 
Top: the experimental average multiplicities are plotted versus the calculated ones. The 
dashed line is the quadrant bisector; well htted points tend to lie on this line. Bottom: 
residual distributions 

Results of hadron multiplicity fit for pp collisions at centre of mass energies 200, 546 and 
900 GeV. Experimental average multiplicities are plotted versus calculated ones. The 
dashed lines are the quadrant bisectors: well fitted points tend to lie on these lines 

Results of hadron multiplicity hts for e’''e“collisions at centre of mass energies 29, 35 
and 44 GeV. Experimental average multiplicities are plotted versus calculated ones. The 
dashed lines are the quadrant bisectors: well htted points tend to lie on these lines 

Results of hadron multiplicity hts for e+e“collisions at centre of mass energy of 91.2 GeV. 
Experimental average multiplicities are plotted versus calculated ones. The dashed lines 
are the quadrant bisectors: well htted points tend to lie on these lines 

Freeze-out temperature values as a function of centre of mass energy. The error bars 
within horizontal ticks for 27.4 GeV pp collisions and for 91.2 GeV e’''e“collisions are 
the ht errors; the overall error bars are the sum in quadrature of the ht error and the 
systematic error related to data set variation (see text) 

Strangeness suppression parameters 7 ^ as a function of centre of mass energy. The error 
bars within horizontal ticks for 27.4 GeV pp collisions and for 91.2 GeV e^e“ collisions 
are the ht errors; the overall error bars are the sum in quadrature of the ht error and the 
systematic error related to data set variation (see text) 

Dependence of htted parameters and primary average multiplicities (top) on the mass 
cut-oh in the hadron mass spectrum for pp collisions at a centre of mass energy of 900 
GeV 
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